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Near the black hole threshold in phase space, the black 
hole mass as a function of the initial data shows the "critical 
scaling" M ~ C(p — f>*) 7 , where p labels a family of initial 
data, p* is the value of p at the threshold, and the critical 
exponent 7 is universal for a given matter model. The black 
hole charge Q obeys a similar law. To complete the picture, 
we include angular momentum as a perturbation. For the 
black hole angular momentum L we find the oscillating be- 
havior L ~ Re [(A + iB)(p - p»Y +iuJ ] . The assumptions of 
the calculation hold for p = p/3 perfect fluid matter, and we 
calculate fj, ~ 0.799 and ui ~ 0.231. 
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a. Introduction There has been much interest re- 
cently in critical phenomena at the black hole thresh- 
old. Take a smooth one-parameter family of smooth, 
asymptotically flat initial data for general relativity (with 
matter) that crosses the black hole threshold, and fine- 
tune the parameter, so as to get close to the black hole 
threshold. For many families, one can make arbitrarily 
small black holes this way. As discovered by Choptuik 
W , the black hole mass near the threshold then scales as 
M ~ C(p — p*) 7 , where p is the parameter of the initial 
data, is its value at the black hole threshold, and 7, 
the "critical exponent" is a universal constant common 
to all one-parameter families, while C depends on the 
family. 

In a nutshell, critical phenomena arise because gen- 
eral relativity, for a given matter model and symmetry, 
admits a "critical solution". By this we understand a 
strong-field solution that is an attractor of codimension 
one sitting on the black hole threshold, and whose basin 
of attraction is the black hole threshold. Perturbatively 
this means that the critical solution has precisely one un- 
stable mode, and that it decays into a black hole, or into 
flat space, depending on the sign of that mode. Criti- 
cal solutions found so far are either static or periodic in 
time (type I), or continuously or discretely self-similar 
(type II). Both can arise in the same matter model 
Here we limit attention to type II, which gives rise to the 
power-law scaling of the mass. For a review of critical 
phenomena, see 

If a small amount of electric charge is present in the ini- 
tial data, the black hole charge Q obeys a similar power 
law as the mass M, with a new critical exponent 8 H||. 
Black holes are completely characterized by their mass 
M, electric charge Q and angular momentum L. To 



complete the picture, here we consider initial data with 
small angular momentum. Our calculation assumes the 
existence, for a given matter model, of a spherically sym- 
metric self-similar solution with exactly one unstable per- 
turbation mode, including nonsphcrical perturbations. 

Recently, we have shown for the first time that a spher- 
ically symmetric critical solution, namely that for perfect 
fluid matter with equation of state p = p/3, has no un- 
stable nonspherical perturbation modes |7]]. The present 
abstract calculation therefore applies to this particular 
matter model, and gives to a quantitative prediction for 
it. 

b. Notation The calculation we give here is hardly 
more than an application of dimensional analysis. There- 
fore we set up a compact notation that hides all technical 
details specific to general relativity, in order to bring out 
the essential scaling argument. Let g M „ be the space- 
time metric in the coordinates (r,x,9, (p), with 9 and tp 
the usual Euler angles, and let surfaces of constant r 
be spacelike. If the conformal metric g^ v = e 2r g^ is 
independent of r, then the spacetime is homothctic, or 
continuously self-similar (CSS). In this case, r plays the 
double role of being the time coordinate and of being the 
negative logarithm of spacetime scale, while x is dimen- 
sionless. A simple example of such a coordinate system 
is t = — ln(— i) and x = r/(— t), where t and r are a 
standard choice of time and radial coordinate. 

Let Z(T,x,9,(p) stand for plus a suitable set of 
conformally rescaled matter variables. (For perfect fluid 
matter, for example, these are p = e~ 2T p for the den- 
sity, and — e T Ufj_ for the four- velocity.) Z(r,x,9,ip) 
then symbolizes a solution of the field equations, while 
Zo(x,9,tp) symbolizes Cauchy data. (To make this no- 
tion precise, we have to transform from g^ v to suit- 
ably rescaled Cauchy data hij and Kij.) Finally, let 
Z(t, x, 9, tp) = Z*(x) stand for the critical solution, which 
is both spherically symmetric and CSS. (We overlook a 
trivial 9 dependence, namely the factor sin 2 9 in g vv , for 
the sake of notation.) 

For a solution that is a general perturbation of a 
self-similar, spherically symmetric background solution 
Z*(x), we can make the ansatz 

00 / 00 

Z(r,x,6,<p)~Z*(x) + J2 E E 

2=0 m=-l n=0 

Re [Ci mn {p) e Ai " T f lmn (x, 9tp)] . (1) 

Because the background is spherically symmetric, we can 
decompose all linear perturbations into spherical har- 
monics labeled by I and m. Because it is independent 
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of r, the perturbation equations of motion separate, and 
the perturbations can be assumed to have a simple expo- 
nential dependence on r, where the separation constants 
Ximn and the mode functions fi mn are complex. We as- 
sume here that the perturbation spectrum is countable 
and discrete, and label it by n for each value of I and m. 
Because Z is real, the Ximn and fimn must be real or form 
complex conjugate pairs. Because of spherical symmetry, 
the Ximn are in fact the same for all m, and are there- 
fore labelled Xi n . The C; mn , also complex, are the free 
parameters of the generic perturbation. Here we want to 
consider one-parameter families of initial data labelled by 
a parameter p. The dependence of the solutions on p is 
then a dependence of the free constants Ci mn on p. The 
fimn correspond to rescaled metric and matter perturba- 
tions in a suitable gauge, such as Regge- Wheeler gauge, 
or in a gauge-invariant framework. We have addressed 
the technical details of gauge-fixing elsewhere . 

c. Review of mass scaling We now consider the time 
evolution of initial data that start close to the spherically 
symmetric self-similar critical solution, so that the expan- 
sion (|]J) applies initially. The critical solution contracts 
to a curvature singularity at r = oo. As the singularity 
is approached from the past, with r increasing, a pertur- 
bation grows (decays) if the real part of its A is positive 
(negative). By definition, a critical solution has precisely 
one unstable mode. Being unique, this unstable mode 
must be spherically symmetric, because all higher values 
of I are (21 + l)-fold degenerate, and its A must be real. 
We label its mode function / o* and its eigenvalue Ao*. 

Close enough to the singularity, we can neglect all 
stable perturbation modes. The unstable perturbation 
keeps growing however, and eventually the solution devi- 
ates nonlinearly from the critical solution. The solution 
eventually forms a black hole for p > p* , but disperses for 
p < p* . As p — + p* from either side, the solution approx- 
imates the critical solution to ever smaller scales, that is 
to ever larger r. This means that Coo*(p) must vanish 
precisely at p — p» . Linearizing Coo* around p = p* , and 
neglecting all decaying perturbations, we have 



Z ~ Z.(x) + ^^(p*)(p -p*)e A °* T /oo*(z)- 
For p > p» , we define t p as a function of p by 
<9Coo* 



dp 



'-(P*)(P - P*)e 



where e is an arbitrary small constant. 

We now consider the Cauchy data Zq given by 

Z(T p ,x,8,ip) ~ Z*(x) + ef ao *(x) = Z (x). 



(2) 



(3) 



(4) 



We have chosen e small enough that the growing mode 
has not yet gone nonlinear (its exact value does not mat- 
ter). For larger r, the linearized ansatz ([!]) breaks down, 
and we know only that the solution forms a black hole 
and settles down to Schwarzschild. We have defined t„ 



so that the rescaled variables Z are independent of p 
in the Cauchy data Z . The field equations themselves 
are scale-invariant (or asymptotically scale-invariant on 
small scales). As a result, there is only one scale, e~ Tp , 
in the problem, and therefore the black hole mass must 
be proportional to it. We have 

M ~ M e~ T ' ~ C(p -p,) 7 . 7 = t^, (5) 

Ao* 

where Mo is a universal constant, and C is a family- 
dependent constant absorbing <9Coo*/dp(p*). This 
derivation of the critical scaling of the black hole mass 
has been given before p|-pl|. Collapse simulations show 
that the critical solution is an attractor, and the mass 
scaling (||) holds, even beyond the strict applicability of 
the perturbation ansatz (gj). 

d. Angular momentum scaling Now we generalize the 
above discussion to initial data with nonvanishing angu- 
lar momentum. Our discussion holds strictly only for 
initial data which are almost spherically symmetric, so 
that deviations from spherical symmetry can be treated 
as linear perturbations throughout the evolution. In par- 
ticular this means that our results hold only for black 
holes whose angular momentum is much smaller than 
their mass, so that the final Kerr solution can be treated 
as a perturbation of Schwarzschild. If the critical so- 
lution turned out to be an attractor beyond the linear 
regime also for non-spherical initial data, our discussion 
and results could then hold also for larger deviations from 
spherical symmetry. 

A Kerr black hole metric with small angular momen- 
tum (L <C M 2 ) can be written as a Schwarzschild metric 
with a linear perturbation proportional to L. In fact, the 
gauge can be chosen so that the shift component of the 
metric is 



9t<p 



L 9 
2- sin0— 
r oO 



0(L 3 /M 6 



(6) 



where Pi is a Legendre polynomial, while all other metric 
coefficients are those of the Schwarzschild metric in the 
usual Schwarzschild coordinates, plus perturbations of 
0(L 2 /M 4 ) | ]l2fl . This result can be verified immediately 
by expanding the Kerr solution in Boyer-Lindquist coor- 
dinates around L = 0. The only perturbation to leading 
order in L therefore has I = 1 angular dependence and 
odd parity in the Regge- Wheeler notation. Here the ori- 
entation of the coordinate axes was chosen so that the 
perturbation is purely m = 0. For a general orientation 
of L relative to the coordinate axes, m = —1,0,1 will 
all appear. As a subdominant effect, we therefore in- 
clude the odd I = 1 perturbations in the data Zq even 
though they are subdominant, with the aim of learning 
something about the behavior of the black hole angular 
momentum near criticality. 

We only keep the dominant odd I = 1 mode, that is 
the one with the largest (least negative) real part of A, 
and denote it by Ai*. Generally, Ai* will not be real, but 
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part of a complex conjugate pair. I = 1 is also threefold 
degenerate, with mode functions /i m * for m — —1,0, 1, 
or the three directions in space. Keeping these modes 
results in a small perturbation of the initial data Zq: 

Z( Tp ,x,e,<p) ~ Z Q {x) + e RcX ^5Z Q (T p ,x,6, V ) : (7) 

where 



l 



5Z = Re ^ C lrn *{ Plf y lm ^- T »h mif {x,e,y) (8) 

m=— 1 

is periodic in r p . As long as its amplitude e ReA i* T p is 
small enough, the perturbation 5Zq evolves as a linear 
perturbation of the spacetime generated by the spheri- 
cally symmetric data Z all way until the solution has 
settled down to Schwarzschild with a small odd-parity 
I = 1 perturbation that turns it into a slowly rotat- 
ing Kerr solution. By the black hole no-hair theorem, 
no other perturbations (except an electric charge of the 
black hole) can survive in the final state. This justifies 
our considering odd I = 1 perturbations but neglecting 
I > 1 and even perturbations altogether. In the linear ap- 
proximation, the black hole angular momentum L must, 
on average in r p , be proportional to e RcA i* T p. Because 
this is a dimcnsionless number, we must have 



L 

M2 



3 ReAi» t p 



(9) 



for the overall scaling of L. But there is an interesting 
modulation due to the fact that SZo depends periodically 
on t p , as we shall see now. 

By linearity, the presence of a perturbation eRe/i mt , 
for some small e, in the initial data Zq must give a rise to a 
small angular momentum component L z /M 2 — eA in the 
final black hole, while e Im/i m * gives rise to L z /M 2 = eB. 
Similar proportionalities apply to Re/i m * and Im/i m * for 
m = ±1, and the x and y components of the black hole 
angular momentum. We absorb the universal, but un- 
known proportionality factors and the family-dependent 
constants Ci m * together into six new constants A and 
B. Putting all the factors together, we obtain the fol- 
lowing final result for the angular momentum vector L 
of the black hole as a function of the initial data, for ini- 
tial data near the black hole threshold and near spherical 
symmetry: 



L = Rc (A + iB)(p-p, 



where 



fj, 



ReAi, 



ImAi* 



— Ao* 



-A, 



(10) 



(11) 



and A and B are six family-dependent constants. Note 
that we always have p, > 2j, because ReAi* < by the 
assumptions of this calculation. 



This is our main result. As p — > p* from above in the 
initial data, the angular momentum vector L of the final 
black hole rotates in space as it decreases. It describes 
a rapidly shrinking ellipse of size (p — p*) M . This behav- 
ior formally resembles a damped isotropic 3-dimensional 
harmonic oscillator, with — ln(p — p*) playing the role of 
"time" , u> the frequency and /i the damping. Of course, 
p — p* is not in any sense a time, but a measure of the 
distance, in phase space, of the initial data set labeled by 
p from the black hole threshold. 

e. The axisymmetric case The result @ is partic- 
ularly surprising if we consider a family of strictly ax- 
isymmetric initial data. Then L can only point along the 
symmetry axis, and the ellipse degenerates into a line: 



L = L z = (p — p^Acosluj ln(p — 



(12) 



where A and c are two family-dependent constants. If the 
family of axisymmetric initial data has angular momen- 
tum that does not vanish at p = p*, then the angular 
momentum of the black hole, as a function of the ini- 
tial data, alternates between parallel and antiparallcl to 
the angular momentum of the initial data. The simple 
power-law scaling L z = (p — p*) M one might naively ex- 
pect would be recovered only if Ai* was purely real. The 
oscillating behavior is sketched in Fig. 1. It is perhaps 
as dramatic an illustration of the "forgetting" of initial 
data in critical collapse as the universality of the critical 
solution itself. 

/. Discrete self- similarity Here we have considered 
the effects of a critical solution that is CSS, both for sim- 
plicity of notation and because the one case that our cal- 
culation is actually known to apply to (the perfect fluid) 
is CSS. It is straightforward, however, to generalize our 
result for the angular momentum scaling to hypothetical 
critical solutions that are DSS instead of CSS. The peri- 
odic dependence of the angular momentum in CSS criti- 
cal collapse is a simple cosine, because it derives from a 
time-dependence expAi*T of the leading I = 1 odd per- 
turbation, with Ai* complex. In DSS critical collapse, 
the background critical solution itself is periodic in r, in 
a nontrivial way, with a period A. Therefore, one would 
expect both kinds of periodicity to arise in L. The calcu- 
lation bearing this out is a straightforward generalization 
of the one for the CSS case, and here we only give its re- 
sult. It is 

L = Rc {A4[ln(p - p*) + c] (A + iB) (p - p,)' I+!W } , 

(13) 

where the only difference is the presence of the complex 
3x3 matrix M.. It is universal and can in principle be 
calculated from the mode functions /i m *. In the CSS 
case it can be assumed to be 1. It is periodic in its 
argument ln(p — p*) with period A/7, c is a family- 
dependent constant. The role of M. is similar to that 
of the fine structure of the black hole mass scaling in 
the DSS case |ll],[R| (where the same parameter c also 
appears). 
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g. Predictions for the p = p/3 perfect fluid We con- 
clude this paper with a quantitative prediction, based 
on numerical work reported elsewhere |Q. The previ- 
ously known critical solution (in spherical symmetry) for 
the perfect fluid with equation of state p = p/3 |la] was 
perturbed, with the result that, at least perturbatively 
around spherical symmetry, it remains the critical so- 
lution in the full phase space. The assumptions of the 
present paper therefore apply to this matter model. The 
values Ao* — 2.785 and Ai* ~ — 0.226 ±0.644i were found 
numerically. The value of Ao* gives rise to a critical ex- 
ponent for the black hole mass 7 = 1/A * ~ 0.359, in 
agreement with the value 0.3558 given previously ||[|. 
Therefore, we predict, for this matter model, for initial 
data near spherical symmetry and near the black hole 
threshold, the oscillatory behavior (|l0|), with /1 ~ 0.799 
and w ~ 0.231. For axisymmetric critical collapse in par- 
ticular (which may be simpler to investigate numeric ally 
than the general case) this reduces to the scaling ( |l2| ) 
with the same fi and uj (as Ai* is not purely real). 
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FIG. 1. Schematic plot of the ADM angular momen- 
tum (thick dashed line) and the angular momentum of the 
black hole (thick full line) in axisymmetric collapse in a 
one-parameter family of axisymmetric data crossing the black 
hole threshold, p is the parameter of the initial data, and p* 
is its value at the black hole threshold. The envelope (thin 
dashed line) is the overall power law. There is an infinite 
number of oscillations as p — » p„ from above. 
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